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ABSTRACT 

We have obtained infrared photometry for the soft x-ray transient GU Mus. We 
present J and K s band light curves modeled with WD98, the newest version of the 
Wilson-Devinney light curve modeling code. Using detailed models for the expected 
ellipsoidal variations due to the non-spherical secondary star, we show that the most 
likely value for the orbital inclination is 54±1.5°. This inclination angle is consistent 
with those previously published, but has a much smaller error. This inclination implies 
a primary black hole mass of 6.95±0.6 M©. While we do not see any evidence for con- 
tamination of our infrared light curves from other sources in the system, a conservative 
model with a contamination level of 15% increases the uncertainty in the inclination 
angle to 54°+f 5 . 

Subject headings: binaries: close — stars: black holes — stars: individual (Nova Muscae 
1991) — stars: low mass — stars: variables: other 



1. Introduction 



Soft X-Ray Transients (SXTs) are binary systems believed to contain a black hole primary 
and low mass secondary star. They display large and sudden x-ray and optical outbursts that are 
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believed to be the result of a sudden, dramatic increase in the mass accretion rate onto the compact 
object. In quiescence, SXTs are very faint at x-ray and optical wavelengths, however, in this state, 
the secondary stars can dominate the system luminosity. 

GU Mus (=Nova Muscae 1991; XN Mus 91; GRS 1124-68) was discovered by the Ginga and 
Granat satellites on January 9, 1991 (Makino 1991; Lund 1991). Since then, it has been studied 
by Remillard et al. (1992), Orosz et al. (1996; hereafter OBMR), Casares et al. (1997; hereafter 
CMCMR), Shahbaz, Naylor, & Charles (1997; hereafter SNC), and references therein. These 
authors have found an orbital period of 10.38 hours, a secondary star spectral type of K4V, and 
a secondary star radial velocity semi-amplitude of 406±7 km s -1 . GU Mus has an implied mass 
function, the minimum mass of the compact primary object, of 3.01±0.15 M (OBMR). This mass 
function combined with the radial velocity curve of the secondary star, established the primary 
object as a dynamical black hole. Both in outburst and quiescence, GU Mus closely resembles the 
prototype black hole SXT, V616 Mon (A0620-00). 

In quiescence, most SXT's light curves reveal ellipsoidal variations ascribed to the secondary 
star. Since the secondary star fills its Roche lobe, the surface area seen by an observer on Earth 
changes as the star orbits the compact object. This changing line-of-sight surface area corresponds 
to a changing apparent brightness. The amplitude of the ellipsoidal variations is determined by the 
orbital inclination angle of the system. By combining the orbital inclination angle with the observed 
mass function, the SXT system parameters can be determined. To this end, we have obtained new 
infrared observations of GU Mus and model its ellipsoidal light curves using a sophisticated light 
curve modeling program. 

Other authors have also attempted to model observed light curves of GU Mus and other SXTs. 
In the case of GU Mus all but one of these papers, SNC, have observed the SXT in the optical. 
If one is searching for the purest ellipsoidal variations, one should observe at a wavelength where 
the secondary star provides the majority of the system luminosity. In the optical, the accretion 
disk and hot spot can contaminate, if not dominate, the system luminosity. Even in quiescence, 
weak and variable accretion might be taking place, and there could be a modest amount of dilution 
of the optical light curve by the accretion disk, hot spot, and possible heating of the secondary 
star (SNC). On the other hand, in SXTs with K and M type companions, the secondary star can 
dominate the quiescent binary's luminosity in the near-infrared. For GU Mus, CMCMR find that 
the secondary star supplies between 85 - 88% of the observed flux in the R band. This percentage 
will be even higher at longer wavelengths. Therefore, observations in the infrared will reveal more 
genuine ellipsoidal variations than observations in the optical. 

The only way to obtain the orbital inclination angle in non-eclipsing binary systems is from 
modeling. Previously published inclination angles for GU Mus range from 39° (lower limit from 
Antokhina & Cherepashchuk 1993) all the way to 83° (upper limit from Hua & Lingenfelter 1993). 
More recently, the inclination has been refined to 54°!^ by SNC, and 60°lg by OBMR. These 
translate into primary mass ranges of 3.8 M — 10.5 M Q , and 5 M Q — 7.5 M Q , respectively. 
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The rather large ranges in the inclination angles and primary masses that have been published for 
GU Mus and other SXT systems, has prompted us to determine whether we could more precisely 
determine the primary mass of the black holes in these SXTs, by using the most recent version of 
the Wilson-Devinney light curve modeling program to model new infrared observations. We have 
observed five SXTs in the infrared, but have chosen to present results on GU Mus first out of our 
sample, since it is the simplest of these systems to model. With no hot spot and very low accretion 
rate, we probably have the best case scenario for modeling relatively uncontaminated ellipsoidal 
variations. 

In section 2, we describe our observations and data reduction process, as well as present our 
infrared photometric light curves. Section 3 describes our infrared light curve modeling procedure 
using WD98, the newest version of the Wilson-Devinney light curve modeling code. We provide 
details of how we chose the relevant input parameters, present the resulting models at J, and K s , 
and discuss the error in the orbital inclination angle. Finally, Section 4 discusses the implications 
of the models, and compares our results to those previously published. 

2. Observations & Data Reduction 

GU Mus was observed using OSIRIS (see Pogge et al. 1999) on the Victor M. Blanco 4- 
m telescope at the Cerro Tololo Inter- American Observatory on 2000 February 20 and 21. On 
February 20, GU Mus was observed from 1:06 to 9:08 UT with the camera at the f/7 plate scale 
(0.161" /pixel). On 21 February from 2:16 to 9:12 UT we observed GU Mus with the camera at 
the f/2.8 plate scale (0.403" /pixel) due to poorer seeing conditions. With an orbital period of 
about 10.38 hours, the first observing session covered two-thirds of an orbital period, while the 
second covered three-quarters of an orbital period. Photometric data were obtained in the OSIRIS 
J (A c =1.215 fim) and K s (A c =2.157 //m) filters. Our observing sequence consisted of a single J 
image at one position, a beam switch, two additional J images, a beam switch back to the original 
position, and then one more J image (an ABBA sequence). We then switched to the K s filter 
and repeated the procedure. Each individual f/7 (f/2.8) J image consisted of 1 frame of 180 (75) 
seconds, while the corresponding f/7 (f/2.8) K s images consisted of 2 (10) coadded frames of 35 
(4) seconds each. 

Dome flats at both plate scales were obtained. Before processing, all of the data were linearized 
using the irlincor package in IRAF with the coefficients supplied in the OSIRIS User's Manual 
(Pogge et al. 1999). After averaging the two images at one position, we subtracted them from the 
average of the two images at the other position. These sky, and bias-subtracted images were then 
flat fielded using the appropriate dome flat. Figure 1 shows a J band image of GU Mus and the 
five nearby field stars used for the differential photometry. 

Aperture photometry was performed on GU Mus and the five nearby field stars shown in 
Figure 1. Using the IRAF phot package, a differential light curve in both J and K s was generated 
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with each f/7 point being the average of the two beam switched images, and each f/2.8 point being 
the average of four beam switched images. Our differential photometric results show that over the 
course of our observations, the comparison stars did not vary more than expected from photon 
statistics. The final J and K s differential light curves of GU Mus, phased to the SNC ephemeris, 
are presented in Figure 2. 

Optical observations of GU Mus were taken with the Cassegrain Focus CCD Imager on the 
0.9 meter telescope at CTIO on 2001 March 14 at 3:40 UT. Data were obtained in the V, R, and I 
bandpasses for the purpose of determining the quiescent optical colors of the system. The 600 second 
exposures were zero corrected and flat fielded before aperture photometry was performed. Landolt 
standard fields were also observed in these optical bands, and were used to correct the apparent 
magnitudes of GU Mus. The resulting colors can be found in Table 1, along with previously 
published optical colors for this sytem. 

3. Modeling the Light Curve of GU Mus 

Before we can begin modeling the light curves of GU Mus, we have to determine a large number 
of parameters. Some of these are intrinsic to the binary system (Sect. 3.1), while others arise from 
the study of stellar atmospheres (Sects. 3.3 & 3.4). 

3.1. System Parameters 

The most important input parameter is the nature of the secondary star. In 1996, OBMR 
compared their observed spectrum of GU Mus to the spectra of various comparison stars, and 
determined a spectral type of K4V±1 for the secondary star. In 1997, CMCMR found that the 
best match to their Ha spectrum was that of a K3V - K4V. Both sets of authors state that the 
secondary star in GU Mus must be slightly evolved since a 'normal' K4V would underfill its Roche 
lobe. The secondary can not be a giant, however, since then it would not fit into an orbit with the 
observed orbital period. 

Knowing the secondary's spectral type does not mean that the mass of the secondary is also 
known. These are most likely not 'normal' stars, as they have survived a supernova explosion and are 
losing mass. Also as just noted, GU Mus probably has a slightly evolved secondary. Determining the 
mass of the secondary is difficult, and provides a large uncertainty to our derivation of the primary 
black hole mass. Even though there is no definitive way to measure the mass of the secondary, 
there are techniques to determine the mass ratio, q, for the system. Orosz et al. (1994) analyzed 
Balmer emission lines to find Ki to use as an indicator of the motion of the primary. They then 
compared Ki to their K2 to obtain a q = K1/K2 = 0.133±0.019. Since their Ki probably arose in 
the accretion disk, it is only a limit to the true radial velocity of the primary. CMCMR calculated 
a system mass ratio of q = 0.128±0.04 from combining their v sini value for the secondary star 
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with geometrical constraints. They were surprised to find that their mass ratio was close to the 
value found by Orosz et al. (1994). Fortunately, the orbital inclination of the system is relatively 
insensitive to the mass ratio (SNC). We will return to this point in section 3.6. 

Interestingly enough, CMCMR did not detect a hot spot in their Ha Doppler map of GU Mus. 
This indicates that the Ha luminosity of the hot spot is smaller than, or equal to, that of the 
accretion disk. CMCMR did detect some emission from the accretion disk; however, and suggest 
that, "Perhaps the mass transfer rate in N. Mus is not sufficient to supply the luminosity of the 
hot-spot." 

3.2. Model Setup 

To model the infrared light curves, we have used WD98, the newest version of the Wilson- 
Devinney light curve program (J. Kallrath 1999, private communication; R. E. Wilson 1999, private 
communication). WD98 is an enhanced version of WD95 (Kallrath et al. 1998), updated with new 
features such as the addition of semi-transparent circumstellar clouds, a simple spectral line profile 
capability for fast-rotating stars, an option to work with either observed times or phases, and the 
conversion of all the variables to double precision. Some of the relevant features of WD98 include: 
Kurucz atmosphere models for numerous wavelengths, a choice of three different limb darkening 
laws which are sensitive to changes in gravity, proximity and eclipse effects, the option for hot or cold 
stellar spots, and several different modes of operation for various system geometries. Descriptions 
of the earlier versions of the program can be found in papers by Wilson and Devinney (1971) and 
by Wilson (1979, 1990, 1993). A recent application of WD95 can be found in Milone et al. (2000). 

Briefly, WD98 works as follows. It takes the photospheres of the stars and divides them up into 
a multitude of surface elements. The amount of light coming from each element is calculated based 
on the binary system input parameters. All of these surface elements are then summed together, 
corrected for the line-of-sight geometry, to create the final light curve. 

There are a number of input parameters needed to generate a model light curve for an SXT such 
as GU Mus. We discuss each of the most important of these parameters in the following subsections. 
We have made use of the best available system parameters. We list the most important wavelength- 
independent input values to WD98 in Table 2, and the wavelength-dependent parameters in Table 
3. Units are shown where appropriate. We ran WD98 in a mode set up to produce a model for a 
semi-detached binary with the secondary star automatically filling its Roche lobe (Mode = 5), and 
the primary having such a large gravitational potential, that it essentially has a radius of zero. We 
now discuss some of the input parameters in detail. 
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3.3. Limb Darkening 

The most important parameter that affects both the shape and the amplitude of the ellipsoidal 
variations is limb darkening. WD98 allows you to choose from several different forms of limb 
darkening: linear, logarithmic, or square-root. The linear law, 

/ A (/x) = I(l)(l-x A (l -//)), 

was first investigated by Milne in 1921. In this equation, I\ is the beam intensity at wavelength A, 
fi is the cosine of the angle between the atmosphere normal and the beam direction, and x\ is the 
limb darkening coefficient. As an alternative to this, Klinglesmith & Sobieski (1970) proposed the 
logarithmic law, 

I x (ji) = J(l)(l - x A (l - n) - yxfJdnfji))), 

where y\ is the non-linear limb-darkening coefficient. Diaz-Cordoves & Gimenez (1992) introduced 
the square-root law, 

l x (n) = J(l)(l - x A (l -fi)- yx(l - y/fl)). 

When compared to ATLAS atmosphere models, the logarithmic law appears to fit UV models the 
best, while the square-root law appears better at infrared wavelengths (Van Hamme 1993). Models 
run by Claret (1998) for very low mass, solar metallicity stars (2000 K < T e ff < 4000 K) indicate 
that the square-root law best describes the intensity distribution in the infrared. We ran test models 
of stars with equal temperature and gravity, and found that the logarithmic and square-root laws 
produced nearly indistinguishable light curves. For the final models presented here, the square-root 
limb darkening law was adopted. 

As shown by Alencar & Vaz (1999) the limb darkening coefficients of stars in close binaries 
can be effected by irradiation. Even though there is no light coming from the primary object in an 
SXT, irradiation by the accretion disk and hot spot might exist. While CMCMR did not detect 
a hot spot in the GU Mus system, they did see a very strong signature from the secondary star. 
They attribute this to a very low mass transfer rate. They also estimated that the secondary star 
contributes 85 - 88% of the system flux in the R band. Since the K4V secondary will dominate 
the system flux to a greater degree in the infrared, we have used the normal, non-irradiated, limb 
darkening coefficients. We present additional evidence for this in section 3.7. 

3.4. Gravity Darkening 

The second most important static parameter that affects the amplitude of the ellipsoidal vari- 
ations is gravity darkening. Gravity darkening (a.k.a. brightening) deals with the localized tem- 
perature of a star's surface. The functional form for gravity darkening, T e ff oc g 13 , is wavelength 
independent, and is not strongly affected by changes in either the mixing length or the composition 
of the star. The amount of gravity darkening depends on how energy is transported through the 
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star, and thus is correlated with the mass of the star. For low mass, convective stars, (3 ps 0.32, 
and for stars with radiative envelopes, (3 ~ 1 (Lucy 1967). For the cool secondary star used in our 
models, we have used a value of (3 = 0.38 according to Claret (2000). He has calculated gravity 
darkening exponents for stars of mass 0.08 M to 40 M Q , and has also found that a more accurate 
method of calculating gravity darkening exponents should include the shape of the star. Until 
that is done, however, we will use the gravity darkening exponent found by Claret (2000) for a 
K4V secondary with T e jj = 4,500 K. This gravity darkening exponent changes very little over the 
estimated spectral type range for the secondary star in GU Mus (K3V — K5V). 

3.5. Other input parameters 

There are a variety of other input parameters with less freedom in their selection. For example, 
we must choose the values for the temperature and monochromatic luminosities of the secondary 
star. For our models, we used a K4V secondary temperature of T e jj=4,500 K, a J band luminosity 
of Lj = 0.311 L Q , and a K s band luminosity of L^ s = 0.421 L Q . 

The atmospheres of cool stars are fairly complicated and the details of their spectral energy 
distributions and any changes in such, as a function of temperature, make for complex modeling 
(cf., Allard et al. 1997). Stellar atmosphere codes have to take into account numerous atomic and 
molecular absorption features which affect the limb darkening coefficients. In order to accurately 
model the limb darkening effects of the secondary star in GU Mus, we have used the Kurucz 
atmosphere models incorporated into WD98. 

Irradiation of the secondary star atmosphere in close binaries such as GU Mus can be very 
important. Harlaftis & Filippenko (2000) used doppler tomography to show evidence of a hot 
spot in QZ Vul (=GS 2000+25) where the accretion stream meets the accretion disk. WD98 can 
calculate reflection/re-radiation of this type of irradiation based on the bolometric albedos of the 
two stars. The expected value for radiative envelopes is unity, while the bolometric albedo for the 
convective secondary star is expected to lie somewhere between 0.5 and 1, based on models run by 
Nordlund and Vaz (Nordlund & Vaz 1990; Vaz & Nordlund 1985). This value is dependent on the 
amount of convection in the star: the smaller the mixing length parameter, a = l/H p , the closer it 
is to a radiative atmosphere, and the higher the bolometric albedo. In the case of GU Mus, we used 
a bolometric albedo of 0.676 based on the average of the albedos given in Table 3 of Nordlund & 
Vaz (1990); however, this parameter was irrelevant since there is little evidence for other significant 
sources of luminosity in the system. 

Our GU Mus J and K s band observations and their corresponding WD98 models for five 
different values of the orbital inclination angle (50°, 53°, 55°, 57°, and 60°), are presented in 
Figures 3 and 4, respectively. The points represent the data while the curves represent the models. 
We present five different inclinations to demonstrate how a change in inclination angle affects the 
resulting light curve in these two band passes. For the models shown in Figs 3 & 4, a \ 2 test gives 
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i = 55° as the best fit to the data. For a fixed set of system parameters, we have found that \ 2 
tests can distinguish between orbital inclinations of ±1°. 

3.6. Sensitivity of the light curve shape to input parameter variation 

We have just shown how we arrived at our input parameters. Many of the parameters discussed 
have simply been set to values consistent with the system parameters of GU Mus. Several of these 
input parameters, however, are either derived from observations of GU Mus (T e ff), or are unknown 
(e.g., the orbital inclination or mass ratio). In order to quantify the sensitivity of the models to 
the variations of these parameters, we ran models with inputs covering a wide range of parameter 
space. 

Varying T e // from 4,400 K (K5V) to 4,800 K (K3V), the observed spectral type range of the 
secondary star, did not produce significant differences in the light curves. Thus for this entire range 
of temperatures, with q = 0.10, the best fit orbital inclination angle is i = 53°. Since we did not 
know the exact mass of either the primary or secondary star in GU Mus, the input values of q 
had to be estimated from previous results. To test the sensitivity of the light curve shape to the 
mass ratio, we ran models with g's covering a very large range: 0.08 < q < 0.13. \ 2 minimization 
tests were performed on all of these models to determine the best fit orbital inclination angle. For 
models with a mass ratio of q = 0.13, we found a best fit inclination angle of i = 55°. For q < 0.1, 
we found that i = 53° gave the smallest x 2 value. Thus, our resulting inclination angle will not be 
very sensitive to q. 

3.7. Sensitivity of the light curve shape to other sources of luminosity in the system 

In the previous sections, we have only considered the light from the secondary star. If there 
are other sources of luminosity in the system, they can affect our result for the orbital inclination 
angle. For example, excess light from a constant source other than the secondary star will dilute 
the ellipsoidal variations. Before we proceed any further, we need to investigate the possibility that 
other sources of luminosity exist in this system. 

As discussed above, CMCMR did not detect a hot spot. They did, however, see evidence for 
emission from a disk. If the accretion disk in GU Mus were to contribute infrared light to the 
system, the true orbital inclination angle would be larger than that derived due to the dilution of 
the light from the secondary star. CMCMR suggested that there was a small (12 - 15%) non-stellar 
contribution to the GU Mus flux in the R band. If we consider a worst case scenario where the 
contamination at J is the same as that observed at R (15%), the derived orbital inclination angle 
is increased by four degrees. Due to the relatively flat nature of model accretion disk spectra (e.g, 
fx oc A~ 2 ; Oke 1977), the contamination will almost certainly be less at J than at R. Assuming this 
15% level of contamination and taking into account the range of mass ratios given in the previous 
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section, we obtain an orbital inclination angle of 54_ 1 5 . 

We now present evidence which suggests that the level of contamination in the GU Mus system 
is significantly smaller than stated above. All of the previous optical studies clearly show unequal 
maxima in the observed light curves of GU Mus (OBMR, for example). This implies that one side 
of the star is heated by the accretion disk, and that contamination is an important factor in the 
modeling of the light curve. The observed infrared light curves presented here have equal maxima, 
and thus we do not see evidence for accretion disk irradiation of the secondary star. This suggests 
a minimal level of heating and/or reflection in the GU Mus system. 

We can also investigate the level of phase-dependent contamination in GU Mus by noting any 
changes in the apparent color of the system throughout its orbit. The infrared color at maximum 
phase will differ from that at minimum phase if either a hotter or cooler source of luminosity exists 
within the system. As can be seen in Figures 3 & 4, the J — K colors agree at both maximum and 
minimum phase, thus again ruling out a significant phase-dependent contamination. 

Given that the expected spectral energy distribution of an SXT at minimum light would be 
more strongly contaminated in the optical portion of the spectrum, we can further examine the 
level of contamination by looking at the optical colors of GU Mus. In Table 1, we present the colors 
of GU Mus and compare the averages to the colors of a K4 dwarf and giant (reddened by Ay = 
0.9 mag). The colors of GU Mus are intermediate to those of the K4 giant and dwarf colors. Thus, 
there is very little evidence for significant blue sources of light in the system. The V — R color is 
slightly redder (~ 10%) than the dwarf or giant, but can be explained by emission from Ha. 

Thus, we do not see significant evidence for emission from the accretion disk itself, or its 
irradiation of the secondary star. But even in the most conservative scenario, the effect on our 
derived orbital inclination angle is < 4°. 



4. Results 

We have explored a wide range of parameter space for GU Mus, and have found that the model 
light curves are relatively insensitive to the input parameters, except for the inclination angle. Even 
changing the mass ratio, q, by a very large amount (i.e., changing the secondary star mass from 
0.56 M to 0.90 M for a primary mass of 6.95 M ), changes the derived orbital inclination angle 
of the system by only two degrees (±1°). We have considered other sources of luminosity in the 
system and even under the most conservative scenario, the derived orbital inclination angle is 
54°1^ 5 . We believe, however, that the infrared light curve is not seriously affected by any such 
contamination. Within the ability of WD98 to realistically model the light curves of SXTs, the 
derived orbital inclination of GU Mus is 5^±_Z.5°. Our derived orbital inclination angle agrees with 
SNC's inclination determined from their H band light curve of GU Mus. 

With the orbital inclination angle determined, we can estimate the mass of the primary. To 
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calculate the mass of the primary object in GU Mus, we need to know three things: the orbital 
inclination, the radial velocity of the secondary star, and the mass of the secondary. The radial 
velocity of the secondary star is K2 = 406±7 km s _1 (OBMR). From the observed range in spectral 
types, the main sequence mass of the secondary star is 0.75±0.05 M (Baraffe et al. 1998). (A 
large uncertainty in finding the mass of the primary remains that associated with the mass of the 
secondary star. As noted by both CMCMR and OBMR, the secondary star in this system may 
be slightly evolved. Even though the mass of the secondary in GU Mus did not have a significant 
effect on our orbital inclination angle determination, it plays a large roll in finding the mass of the 
primary.) Using these values, we calculate a black hole mass of 6.95±0.6 Mq. This mass is within 
the range reported by both CMCMR and SNC. 

For an average primary mass of 6.95 Mq, and an average secondary mass of 0.75 Mq, we obtain 
a q = 0.108. Using this mass ratio and a semi-major axis of 4.6 Rq, we find the resulting Roche 
lobe radius of the secondary star in GU Mus to be 0.97 R Q . Note that a main sequence K4V would 
have a radius of 0.73 R , thus filling ~ 75% of the derived Roche lobe radius. Interestingly, both 
the derived radius and optical colors are consistent with the secondary star being a subgiant. Using 
R = 0.97 Rq and T e jj = 4,500 K (appropriate for a K4V), we derive a bolometric luminosity of 
L_Bo2 = 0.34 L Q . Using our infrared photometry, corrected for reddening, we determine a distance 
to GU Mus of 5.1 kpc. 

The development of a consistent predictive theory of accretion flows around black holes is 
one of the main goals of high-energy astrophysics. The origin of the outbursts of SXTs can be 
explained by the accretion disk limit cycle mechanism that was developed to describe the outbursts 
of cataclysmic variables (Cannizzo 1998). Modeling of the outbursts of SXTs, from the x-ray 
through optical, has been attempted by a number of researchers. The consensus appears to be 
that the optical and soft x-ray luminosities are generated within an accretion disk. The hard x-ray 
spectra that develop later in the outburst have been explained as advection dominated accretion 
onto the black hole. Depending on the model chosen to explain the multi-wavelength luminosities, 
the inclination angle can play an important role in the appearance of the outburst. 

For example, Esin, McClintock, and Narayan (1997) calculated theoretical x-ray models to fit 
to various phases that GU Mus exhibited during its outburst. They combined a normal accretion 
disk with an advection dominated accretion flow (ADAF). In their models, they input a black hole 
mass of 6 M , an inclination angle of i = 60°, and a distance of 5 kpc, taken from OBMR, CMCMR, 
and SNC. They found that the during the outburst of GU Mus, both the accretion disk and ADAF 
remained optically thin. If so, the inclination angle does not play a serious role in what is observed. 
But as Esin et al. (2000) note, even if the ADAF is optically thin, the emission will preferentially 
escape along the polar axes, and thus the inclination angle may affect the observed x-ray flux. In 
contrast, Misra (1999), uses a similar model (with identical values of Mi, i, and distance), but one 
where both the ADAF and the accretion disk are optically thick. If this model is correct, then the 
inclination angle plays an important role in what is observed. New models that are constructed to 
explain the outburst of GU Mus should take into account our new results for the mass of the black 
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hole, the orbital inclination angle and the distance. 

5. Conclusions 

We have obtained new J and K s band light curves covering a complete orbit for GU Mus. 
We have used WD98 to model the resulting light curves, and obtained an orbital inclination of 
54±1.5°. This leads to a primary mass estimate of 6.95±0.6 Mq. In quiescence, GU Mus is a 
relatively simple system to model: There is no hot spot, and very little accretion to contaminate 
its infrared light curves. There is also a consensus on the orbital inclination angle for GU Mus, 
unlike the other SXTs. Thus, it provides the opportunity to test our modeling procedure before 
attempting more complex systems. 

Before we can fully understand the outbursts of SXTs, we need to determine values for the most 
important systemic parameters. Unfortunately, it is much harder to determine system parameters 
for non-eclipsing binaries, than for their eclipsing counterparts. The only way to determine the 
orbital inclination in a non-eclipsing system is to model its ellipsoidal variations. WD98 allows 
us to model these variations using the most up-to-date software package available. By modeling 
infrared data, where the contamination from non-stellar sources is minimized, we have a robust 
technique to determine important orbital parameters. 
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Fig. 1.— OSIRIS J band image of GU Mus (a 2 ooo = 11:26:26.7, <5 2 ooo = -68:40:32.6) taken on 
February 20, 2000 with the Victor M. Blanco 4 meter telescope at Cerro Tololo Inter-American 
Observatory. The field is 1.5' x 1.5' with a scale of 0.16" /pixel. The SXT is labeled with a "G" . 
The five comparison stars used for the differential photometry are numbered. 
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Fig. 2. — GU Mus J band (top panel) and K s band (bottom panel) light curves. The data were 
obtained on February 20 (triangles) and February 21 (circles), with OSIRIS on the CTIO Victor M. 
Blanco 4-m telescope. The data are plotted over two phase cycles for clarity. Here and throughout 
this paper we phase our heliocentric corrected data to the Shahbaz et al. (1997) ephemeris. 
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Fig. 3. — J band light curve with five WD98 models. The data are the same as in Figure 2. The 
lines are WD98 models of differing orbital inclination angles: 50° (long-dashed), 53° (solid), 55° 
(dotted), 57° (short-dashed), and 60° (dot dashed). With q = 0.13, a \ 2 minimization finds the 
best fit model to be 55°. 
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Fig. 4.— 



K s band light curve with the five WD98 models as in Figure 3. 
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Table 1. Quiescent Infrared and Optical Colors of GU Mus 



Reference 



V 



B - V 



V - R 



V - I 



J - K 



Year a 



20.5±0.1 
20.66±0.03 
20.55±0.05 
20.35±0.03 
20.83±0.06 



1.60 



1.00 
0.91 



0.90±0.07 



1.66 
1.53 
1.55 
1.83±0.07 



18.14±0.02 1.06±0.10 



1992 
1992 
1992 
1993 
2000/2001 



Mean Observed Color 
K4V (A v = 0.9 mag) 
K4III (A v = 0.9 mag) 



1.60±0.06 
1.36 
1.70 



0.93±0.05 
0.77 
0.88 



1.64±0.14 
1.53 
1.87 



1.06±0.10 
0.83 
1.03 



J Year the data were taken 



References. — (1) Delia Valle et al. 1998; (2) King et al. 1996; (3) Remillard et al. 1992; (4) OBMR; (5) This 
Paper 



Table 2. Wavelength Independent WD98 Input Parameters for GU Mus 



Parameter Value 



Orbital Period (days) 0.432604 

Ephemeris a (HJD phase 0.0) 2448812.669 

Semi-Major Axis (R ) 4.6 

Orbital Eccentricity 0.0 

Temperature of K4V Secondary (K) 4,500 

Mass Ratio (M 2 /Mi) 0.108 

Atmosphere Model Kurucz 

Limb Darkening Law Square-root 
Secondary Star Gravity Darkening Exponent ,9=0.38 

Secondary Star Bolometric Albedo 0.676 



a From Shahbaz, Naylor, and Charles (1997) 
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Table 3. Wavelength Dependent WD98 Input Parameters for GU Mus 



Parameter 


J 




Secondary Star Monochromatic Luminosity (Lq) 


0.311 


0.421 


Secondary Star Square-root Limb Darkening Coefficients: 








0.110 


-0.116 


y\ 


0.531 


0.724 



